Resonant Coherent Phonon Spectroscopy of Single- Walled Carbon Nanotubes 
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Using femtosecond pump-probe spectroscopy with pulse shaping techniques, one can generate and 
detect coherent phonons in chirality-specific semiconducting single-walled carbon nanotubes. The 
signals are resonantly enhanced when the pump photon energy coincides with an interband exciton 
resonance, and analysis of such data provides a wealth of information on the chirality-dependence of 
light absorption, phonon generation, and phonon-induced band structure modulations. To explain 
our experimental results, we have developed a microscopic theory for the generation and detection 
of coherent phonons in single- walled carbon nanotubes using a tight-binding model for the electronic 
states and a valence force field model for the phonons. We find that the coherent phonon amplitudes 
satisfy a driven oscillator equation with the driving term depending on photoexcited carrier density. 
We compared our theoretical results with experimental results on mod 2 nanotubes and found that 
our model provides satisfactory overall trends in the relative strengths of the coherent phonon signal 
both within and between difi'erent mod 2 families. We also find that the coherent phonon intensities 
are considerably weaker in mod 1 nanotubes in comparison with mod 2 nanotubes, which is also in 
excellent agreement with experiment. 

PACS numbers: 78.67.Ch, 63.22.-fm, 73.22.-f, 78.67.-n 



I. INTRODUCTION 

Single-walled carbon nanotubes (SWNT) can be 
viewed as rolled up sheets of graphene, having a 
one-dimensional band structure with unique electronic, 
mechanical, and optical properties. Their electronic 
properties vary significantly, depending on their chi- 
rality indices (n,m), and can be either metallic or 
semiconductingiii^i^i^ Although there are currently 
world-wide efforts to achieve single-chirality samples, a 
standard for fabrication of such samples has yet to be 
established. Resonant Raman spectroscopy (RRS) or 
photoluminescence excitation spectroscopy (PLE) is usu- 
ally used to study chirality-dependent electronic and vi- 
brational properties. However, carbon nanotube sam- 
ples typically contain ensembles of nanotubes with dif- 
ferent chiralities, and the unknown relative abundances 
of different-chirality tubes in such samples often makes 
it challenging to extract reliable parameters on chirality- 
dependent properties from experimental results. 

Resonant Raman spectroscopy can be used to study 
chirality-dependent electron-phonon coupling in nan- 
otubes and can be used to uniquely determine the chi- 
rality of individual tubesi^i^iSiiSiiiS Raman spectroscopy is 
a sensitive probe of ground-state vibrations but is less 
suitable for studying excited state vibrational proper- 
ties. Recently, excited state lattice vibrations in car- 
bon nanotubes have been studied with coherent phonon 
(CP) spectroscopy .^^i^^'^'^i^^ In CP spectroscopy, coher- 
ent phonon oscillations are excited by pumping with 
an ultrafast pump pulse and are detected by measur- 



ing changes in the differential transmission using a de- 
layed probe pulse. The CP intensity is then obtained by 
taking the temporal power spectrum of the differential 
transmission. The peaks in the power spectrum corre- 
spond to coherent phonon frequencies. Coherent phonon 
spectroscopy allows direct measurement of excited state 
phonon dynamics in the time domain including phase in- 
formation and dephasing times. 

We have developed a technique that allows us to 
study chirality dependent properties of nanotubes in an 
ensemble.— This is described in Section [TTl By shaping 
the pump pulse, we incorporate quantum control tech- 
niques in CP spectroscopy. Using pre-designed trains 
of femtosecond optical pulses, we have selectively ex- 
cited and probed coherent lattice vibrations of the radial 
breathing mode (RBM) of specific chirality single- walled 
carbon nanotubes. We are able to gain information on 
light absorption, coherent phonon generation, and co- 
herent phonon-induced band structure modulations. We 
find that coherent RBM phonons can be selectively ex- 
cited by using a train of pump pulses whose repetition 
rate is in resonance with the desired phonon frequency. 
By exciting only those phonon modes with a specific fre- 
quency, we can selectively study nanotubes with the same 
chirality in an ensemble of tubes. 

In order to explain our experimental results, we de- 
velop in Section IIIII a microscopic theory for the genera- 
tion and detection of coherent phonon lattice vibrations 
in carbon nanotubes by ultrafast laser pulses. We use 
a third nearest neighbor extended tight-binding (ETB) 
model to describe the electronic states over the entire 
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Brillouin zone while the phonons are treated in a valence 
force field model. In treating the electrons and phonons, 
we exploit the screw symmetry of the nanotube to drasti- 
cally simplify the problem. Equations of motion for each 
CP vibrational mode arc obtained, using a microscopic 
description of the electron-phonon interaction based on 
direct evaluation of the three-center electron-phonon ma- 
trix elements using ab initio wavefunctions and screened 
atomic potentials. For each CP active mode we find that 
the CP amplitudes satisfy a driven oscillator equation 
with a coherent phonon driving function that depends 
on photoexcited hot carrier distributions. An ultrafast 
laser pulse generates electron-hole pairs and the driving 
function rises sharply in a step-like fashion. If the pulse 
duration is shorter than the phonon oscillation period, 
the rapid initial jump in the coherent phonon driving 
function gives rise to oscillating coherent phonon ampli- 
tudes. 

Carbon nanotubes with the same values of 2n -I- m are 
said to belong to the same family (with index 2n + to). 
Carbon nanotubes in a given family are metallic if 
mod (n — TO, 3) = and semiconducting otherwise. The 
semiconducting tubes are classified as either mod 1 or 
mod 2 depending on whether the value of mod (n— to, 3) 
is 1 or 2. In CP spectroscopy, we find that a strong signal 
is obtained when we pump at the allowed nanotube En 
optical transitions. We found experimentally that, for 
the RBM modes, the CP intensity within a mod 2 family 
tends to decrease with chiral angle and the decrease in 
CP intensity with chiral angle is found to be much more 
pronounced for the En feature. We also found that CP 
intensities are considerably weaker in mod 1 families in 
comparison with mod 2 families. In general, the E22 
CP intensities in mod 2 families are stronger than the 
Ell features while the opposite is true in mod 1 tubes. 
For RBM modes in mod 1 tubes, the En CP intensities 
tend to decrease with increasing chiral angle within a 
given family. As the family index (2n-|-TO) increases, the 
Ell CP intensity in mod 1 tubes decreases. Finally, we 
compared our theoretical results with experimental CP 
spectra in mod 2 nanotubes and found that our theoreti- 
cal model correctly predicts the experimentally observed 
overall trends in the relative strengths of the CP signal 
both within and between mod 2 families. We found dis- 
crepancies between our theoretical predictions with re- 
gard to the peak positions and lineshapes. These dis- 
crepancies can be qualitatively attributed to Coulomb 
interactions which have not yet been included in the cal- 
culations. 

We do not consider the Coulomb interaction and exci- 
tonic effects in our theoretical model for reasons of sim- 
plicity and tractability. It has been pointed out that 
excitonic effects are important for understanding the op- 
tical properties of small diameter carbon nanotubesj^ 
In a number of Raman scattering theories the Coulomb 
interaction is neglected, but nevertheless the computed 
RBM Raman spectra can explain many experimental 
measurementsJiiSiii^ It is worth dwelling a little on the 



reason for this. Jiang et al<^ recently undertook a study 
of the exciton-photon and exciton-phonon matrix ele- 
ments in single-walled carbon nanotubes using a tight- 
binding model. These authors found that for the RBM 
and G-band modes, the phonon matrix elements in the 
exciton and free particle pictures are nearly the same. 
However, values for the exciton-photon matrix elements 
are on the order of 100 times greater than the electron- 
photon matrix elements computed in the free particle 
picture. Thus, when we discuss the photoexcitation of 
carriers, the actual photoexcited carrier densities will be 
different from what we predict. On the other hand when 
we discuss the dependence of the coherent phonon ampli- 
tudes on the Eii transition energies or the tube chirail- 
ity, the present discussion has a physical meaning. Thus 
we expect reasonable agreement between experiment and 
theory for relative coherent phonon amplitudes and the 
relative strengths of the computed CP spectra. However, 
the peak positions and line shapes of the CP spectra will 
be altered by the neglect of excitonic effects. 

Apart from being useful in resonant CP spectroscopy 
of chirality specific nanotubes, we note that laser induced 
coherent phonons in carbon nanotubes may also have im- 
portant practical applications in the fabrication of car- 
bon nanotube electronic devices. In recent years laser in- 
duced coherent phonon generation has been theoretically 
studied using molecular dynamics techniques j^i^ii^iS^i^ 
Garcia et al^ have simulated laser induced coherent 
phonons in a mod 1 zigzag (10,0) capped carbon nan- 
otube using a formalism that combines a nonadiabatic 
molecular dynamics method and a density matrix ap- 
proach to describe the dynamics of the carbon ions and 
valence electrons. DumitricS et ali^^^ have theoretically 
studied the possibility of achieving selective cap opening 
in (10,0), (5,5) and (8,4) capped carbon nanotubes driven 
by laser induced coherent phonons using nonadiabatic 
molecular dynamics simulations based on a microscopic 
electronic model. 

It is well known that self-assembled carbon nanotubes 
suffer from structural imperfections that modify their 
electronic, optical, and mechanical properties. Such de- 
fects pose a problem for the fabrication of nanotube 
based electronic devicesi^ A common type of defect in 
carbon nanotubes is the (5-7) pair defect introduced 
by applying a Stone- Wales transformation to the nan- 
otube structure;^ In the Stone- Wales transformation, 
four hexagons are replaced by two pentagons and two 
heptagonsi^ The possibility of eliminating such defects 
using laser generated coherent phonons has been stud- 
ied theoretically in armchair and zigzag nanotubes using 
nonabiabatic molecular dynamics simulations. Romero 
et al?"^ have studied the response of armchair nanotubes 
with (5-7) pair defects to ultrafast laser pulses and found 
that when the fraction of photoexcited electrons exceeds 
a critical threshold ( around 7% ) the resulting coherent 
phonon oscillations cause the nanotube to undergo an in- 
verse Stone- Wales structural transition that heals the de- 
fect. More recently Valencia et al?^ have extended these 
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FIG. 1: (color online) Experimental setup for chirality- 
selective excitation of coherent lattice vibrations in single- 
walled carbon nanotubes through ultrafast pulse shaping. 
The inset shows an example of multiple pulse trains with a 
repetition rate of 6.15 THz and a central wavelength of 800 nm 
(1.55 eV). 



studies to zigzag nanotubes and found similar results. 

Recently, Jeschke et al^ have theoretically studied the 
structural response of nanotubes of different chiralities 
to femtosecond laser excitation using molecular dynam- 
ics simulations. They found that carbon nanotubes may 
transform into more stable structures under the appropri- 
ate conditions. Such investigations may be important for 
technological applications. For example, nanotubes ex- 
cited by lasers above a certain threshold may tear open 
and interact with other tubes leading to the creation of 
new structures. 



II. EXPERIMENT 

In this section, we demonstrate the implementation 
of ultrafast pulse-shaping to excite the coherent radial 
breathing modes of specific chiralities, providing a defini- 
tive ability to study single chirality nanotubes from an 
ensemble sample. Our method exploits selective excita- 
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FIG. 2: (color online) Generation and detection of coherent 
phonons of the radial breathing mode in single- walled carbon 
nanotubes. (a) Time-domain transmission modulations due 
to coherent RBM vibrations in ensemble single-walled carbon 
nanotube solution that were generated by standard pump- 
probe spectroscopy without pulse shaping, (b) Fourier trans- 
formation of time-domain oscillations with chirality assigned 
peaks. 



tion to not only extract chiral-dependent band gap mod- 
ulations, but also utilizes information from the probe 
energy dependence of the phase and the amplitude of 
the coherent phonon oscillations to reconstruct excitation 
profiles for the E22 transitions. In particular, our ob- 
servation of probe-energy-dependent phase reversal pro- 
vides direct, time-domain evidence that for coherent ra- 
dial breathing modes the band gap oscillates in response 
to the nanotube diameter oscillations. 

The sample used in this study was a micelle-suspended 
SWNT solution, where the single-walled carbon nan- 
otubes (HiPco batch HPR 104) were suspended as in- 
dividuals with sodium chelate. The optical setup was 
that of standard degenerate pump-probe spectroscopy, 
but chirality selectivity of RBM oscillations was achieved 
by using multiple pulse trains, with a pulse-to-pulse in- 
terval corresponding to the period of a specific RBM 
mode.^^ Among different species of nanotubes, those hav- 
ing RBM frequencies that are matched to the repetition 
rate of multiple pulse trains will generate large amplitude 
coherent oscillations with increasing oscillatory response 
to each pulse, while others will have diminished coherent 
responses Jii^iii^ The tailoring of multiple pulse trains 
from femtosecond pulses was achieved using the pulse- 
shaping technique developed by Weiner and Leaird.^^ As 
depicted in Fig. [U pulse trains are incident on an en- 
semble of nanotubes as a pump beam, whereas coherent 
RBM oscillations are monitored by an unshaped, Gaus- 
sian probe beam. 

Real-time observation of coherent RBM oscillations 
is possible without pulse-shaping by employing stan- 
dard femtosecond pump-probe spectroscopyi ^^i^^'^'^ Fig- 
ure [Ha) shows transmission modulations of the probe 
beam induced by coherent lattice modulations, which 
were generated by pump pulses with a pulse width of 50 fs 
and a central wavelength of 800 nm (1.55 eV). The time- 
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FIG. 3: (color online) (Left) Time- domain coherent RBM 
oscillations selectively excited by multiple pulse trains via 
pulse-shaping with corresponding repetition rates from 6.15- 
7.07 THz. (Right) Fourier transformations of corresponding 
oscillations, with their dominant nanotube chirality (n,m) in- 
dicated. 



domain beating profiles reflect the simultaneous genera- 
tion of several RBM frequencies from nanotubes in the 
ensemble with different chirahties, which are clearly seen 
in Fig. m^b) with the Fourier-transformation of the time- 
domain data. Although resonance conditions and mode 
frequencies lead to the assignment of chiralities to their 
corresponding peaksfi^ obtaining detailed information on 
dynamical quantities such as the phase information of 
phonon oscillations becomes rather challenging. Addi- 
tionally, if adjacent phonon modes overlap in the spectral 
domain, this can lead to peak distortions. 

However, by introducing pulse-shaping, multiple pulses 
with different repetition rates are used to excite RBM 
oscillations, and as shown in Figs. [3I^a)-[3Jd), chirality se- 
lectivity was successfully obtained. With the appropriate 
repetition rate of the pulse trains, a single, specific chiral- 
ity dominantly contributes to the signal, while other nan- 
otubes are suppressed. For example, by choosing a pump 
repetition rate of 7.07 THz, we can selectively excite only 
the (11, 3) nanotubes, as seen in Fig.[3](a). Similarly, with 
a pump repetition rate of 6.69 THz, the (10, 5) nanotubes 
are selectively excited, as seen in Fig.[3l^b). The accuracy 
of selectivity depends on the number of pulses in the tai- 
lored pulse train as well as the distribution of chiralities 
in the nanotubes ensemble. Furthermore, selective exci- 
tation of a specific chirality also requires the pump energy 
to be resonant with the corresponding £'22 transition for 
each chirality-specific nanotubes. 
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FIG. 4: (color online) Differential transmission as a function 
of time delay at probe wavelengths of 780 nm, 795 nm, and 
810 nm for the selective RBM excitation of the (11,3) nan- 
otubes. There is a tt phase shift between the 780 nm and 
810 nm data. These three wavelengths (from the top to the 
bottom of the figure) correspond to photon energies above, 
at and below the energy of the second exciton resonance, re- 
spectively, of (11,3) nanotubes. 



The ability to excite single-chirality nanotubes allows 
us to perform detailed studies of excited states of single- 
walled carbon nanotubes. For example, by placing a se- 
ries of 10-nm band pass filters in the probe path before 
the detector, we can measure the wavelength-dependence 
of RBM-induced transmission changes in order to under- 
stand exactly how the tube diameter changes during co- 
herent phonon RBM oscillations and how the diameter 
change modifies the nanotube band structure. As seen 
in Fig. m the differential transmission is shown for three 
cases, from top to bottom, corresponding to probe pho- 
ton energies above-resonance, on-resonance, and below- 
resonance, respectively, for selectively-excited (11,3) car- 
bon nanotubes. Although the transmission is strongly 
modulated at the RBM frequency (7.07 THz) for all 
three cases, the amplitude and phase of oscillations vary 
noticeably for varying probe wavelengths. Specifically, 
the amplitude of oscillations becomes minimal at reso- 
nance, and, in addition, there is clearly a 7r-phase shift 
between the above- and below-resonance traces. Because 
the band gap energy and diameter are inversely related 
to each other, and because it is the RBM frequency at 
which the diameter is oscillating, we can conclude from 
this data that the energy of the E22 resonance is oscil- 
lating at the RBM frequency. Namely, when the band 
gap is decreasing, absorption above (below) resonance is 
decreasing (increasing), resulting in positive (negative) 
differential transmission. 

We can also look at the short response to see how the 
diameter changes in response to ultrafast excitation of 
electron- hole pairs by the pump pulse. In Fig. OJa), we 
plot the differential transmission data taken at 780 nm 
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FIG. 5: (color online) (a) Differential transmission as a func- 
tion of time delay together with the pump pulse train (dark 
blue or dark gray) as well as fit (light blue or light gray up to 
4ps) to exponentially-decaying sinusoidal oscillations. (b,c) 
Data near time zero for two wavelengths corresponding to 
energies above and below the second exciton resonance, re- 
spectively, of (11,3) nanotubes. 



(1.59 eV) together with the pump pulse train, with time 
zero corresponding to the center of the pulse train. Here 
we note that an increase in the absorption corresponds to 
a decrease in the differential transmission. Figures O^b) 
and [3c) show data near time zero for two wavelengths 
corresponding to energies above and below the second ex- 
citon resonance, respectively, of (11,3) nanotubes. The 
sign of the differential transmission oscillations in the first 
quarter-period, where the time delay varies from 0.0 ps 
to 0.07 ps, is positive (negative) for the above (below) 
resonance probe, indicating that there is an initial de- 
crease (increase) in absorption for energies above (below) 
resonance, demonstrating that the diameter of the nan- 
otube initially expands, taking into account the fact that 
the resonance energy is inversely related to the diameter. 
This initial expansion of the tube diameter is in agree- 
ment with our theoretical predictions for the photoexci- 
tation of coherent phonon REM oscillations by ultrafast 
laser pulses pumping near the E22 transition energy in 
mod 2 nanotubes [e.g., (11,3) tubes]. 



III. THEORY 

We have developed a microscopic theory for the gen- 
eration of coherent phonons in single-walled carbon nan- 
otubes and their detection by means of coherent phonon 
spectroscopy experiments. Our approach is based on ob- 



taining equations of motion for the coherent phonon am- 
plitudes from the Heisenberg equations of motion as de- 
scribed by Kuznetsov et al. in Ref. [l^. In our theoretical 
model, we explicitly incorporate the electronic energies 
and wavefunctions for the tt electrons, the phonon dis- 
persion relations and the corresponding phonon modes, 
the electron-phonon interaction, the optical matrix el- 
ements, and the interaction of carriers with a classical 
ultrafast laser pulse. For simplicity, and to make the 
problem tractable, we neglect the many-body Coulomb 
interaction and interactions with the surrounding liquid 
medium in the micelle-suspended nanotube ensemble. 

We are able to treat nanotubes of arbitrary chirality 
by exploiting all the screw symmetry operations. This 
allows us to examine trends in the CP signal strength 
within and between nanotube families. In addition, we 
gain something in our conceptional understanding by de- 
riving a simple driven oscillator equation for the coher- 
ent phonon amplitudes where the driving function de- 
pends explicitly on the time-dependent photoexcited car- 
rier distribution functions. In the limit where we ignore 
Coulomb interactions, the driven oscillator equation for 
the coherent phonon amplitudes turns out to be exact 



A. Electron Hamiltonian 

We treat carbon nanotube tt and tt* electronic states in 
the extended tight-binding (ETB) formalism of Porezag 
et aZ.^- In the ETB model, the tight-binding Hamilto- 
nian and overlap matrix elements between tt orbitals on 
different carbon atoms are functions of the interatomic 
distance. The position-dependent Hamiltonian and over- 
lap matrix elements are obtained from a parametrization 
of density-functional (DFT) results in the local-density 
approximation (LDA) using a local orbital basis set as 
described in Ref. 29. Our computed energy dispersion 
relations for the bonding tt and anti-bonding tt* bands in 
graphene are plotted in Fig. [S] along high symmetry di- 
rections in the hexagonal two-dimensional Brillouin zone. 
For comparison, we also plot the graphene energy dis- 
persion relations obtained from the simple tight-binding 
model (STB), described in Ref. H, in which only nearest 
neighbor Hamiltonian and overlap matrix elements are 
considered. In the STB model, the values of the nearest 
neighbor Hamiltonian and overlap matrix elements are 
-3.033 eV and 0.129, respectively^. The ETB and STB 
models agree with each other near the K and K' points 
in the Brillouin zone and in carbon nanotubes these are 
the states that give rise to the low lying conduction and 
valence subbands that we are interested in. 

In a carbon nanotube with chiral indices (n,m), a 
translational unit cell can be found such that the atomic 
structure repeats itself after translation of the tube by a 
translational vector T parallel to the tube axist^ The re- 
sulting Brillouin zone is one-dimensional with |fc| < tt/T 
and the number of two-atom hexagonal cells in each 



6 




.10-1 1 1 1 

r M K r 



FIG. 6; (color online) Extended tight-binding energy disper- 
sion relations for the bonding n and anti-bonding tt* bands in 
graphene along high symmetry directions in the Brilloin zone 
are shown as solid black lines. Energy dispersion relations in 
the nearest neighbor tight-binding model described in Ref. Tl 
are shown as dash-doted red lines. 



translational unit cell ia* 

^ 2{n'^ + nm + m?) 

hex - gcd(2n + m, 2m + n) 

where gcd(i, j) is the greatest common divisor of integers 
i and j. If we only make use of translational symmetry in 
formulating the electronic problem, the resulting size of 
the Hamiltonian and overlap matrices is 2A^j^g^ x iN-^^^^ 
if we retain one tt orbital per site. 

In practice, the size of the electronic Hamiltonian and 
overlap matrices obtained using the nanotube transla- 
tional unit cell can become prohibitively large, especially 
for chiral nanotubes with n > m ^ 0. Fortunately, we 
can reduce the size of the electronic problem by further 
exploiting the symmetry of the nanotube. As pointed out 
in Ref. [iTl a two-atom hexagonal unit cell in graphene, 
with the two carbon atoms labeled A and B, can be 
mapped onto the nanotube by applying two different 
screw operations. If we make use of the screw symmetry 
operations, we can block diagonalize the 2N^^^y^ x 2Ny^^^ 
Hamiltonian and overlap matrices into 2x2 subblocks 
which we label /i. In carbon nanotubes, the subblock in- 
dex /J, labels the cutting lines in the zone folding picture. 
For states near the Fermi energy, the cutting line num- 
bers ^ have a nice geometrical interpretation as pointed 
out in Ref. [s^. 

Derivations for the Hamiltonian and overlap matrices 
for the electronic states in a carbon nanotube are given 
in Appendix[Xl We let s = v,c label the valence and con- 
duction band states, and the electronic energies Es^{k) 
for a given cutting line are obtained by solving the matrix 
eigenvalue problem in Eq. (|A8|) . The second quantized 



electron Hamiltonian is simply 

where c^^^ creates an electron in the state with energy 

B. Phonon Hamiltonian 

Following Jiang et al^ and Lobo et al^, we treat lat- 
tice dynamics in a carbon nanotube using a valence force 
field model. In our force field model, we include bond 
stretching, in-plane bond bending, out-of-plane bond 
bending, and bond twisting potentials. In constructing 
the valence force field potentials, we take care that they 
satisfy the force constant sum rule which requires that 
the force field potential energy remain invariant under 
rigid translations and rotations (see Ref. [H, p. 131). As 
pointed out in Ref. IsS , a number of calculations in the 
literature use force field models that violate the force con- 
stant sum rule and, as a result, fail to reproduce the long 
wavelength flexure modes predicted by elasticity theory. 

Our valence force field model as described in Ap- 
pendix |B] has seven force constants, four due to bond 
stretching interactions out to fourth nearest neighbor 
shells and one each from the remaining three interac- 
tions. To determine these seven force constants, we fit 
our model results for planar graphene to the model of 
Jishi et al.^^ Our best fit dispersion relations are shown 
in Fig. [7] as solid lines while the results of Jishi et al. are 
shown as solid circles. 

We should point out that since our force field model 
contains force constants that are independent of the den- 
sity of photoexcited carriers, it cannot describe phonon 
softening which is observed at high values of the laser 
fluence. However, a rather high value of the laser flu- 
ence is generally needed to generate a high density of 
photoexcited carriers. In the case of metallic nanotubes, 
the chirality dependent frequency shift of the RBM and 
G modes has been studied by Sasaki et al. in Refs. [H 
and as a function of the Fermi energy. In the case 
of the RBM mode they find that armchair nanotubes do 
not exhibit any frequency shift while zigzag nanotubes 
exhibit phonon softening)^ 

The phonon energies and corresponding mode displace- 
ment vectors are obtained by diagonalizing the dynam- 
ical matrix. In graphene, there are two atoms per unit 
cell giving rise to the six phonon modes shown in Fig. [T] 
In a carbon nanotube, the size of the dynamical matrix 
is SA^jj^g^ X 6A^j^g^. By making use of the nanoutube 
screw symmetry operations, we can block diagonalize the 
dynamical matrix into 6x6 subblocks which we label 
V = . . . Af^gj^ — 1. Again, the subblock index v labels 
the cutting fines. In a carbon nanotube, the phonon en- 
ergies hujf}i,{q) are obtained from solving the dynamical 
matrix eigenvalue problem in Eq. (jBlOp where v is the 
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FIG. 7: Graphene phonon dispersion relation. The solid 
circles are the valence force field model of Jishi et al.'^^ and 
the solid lines are the best fit results for our seven parameter 
valence force field model. 

cutting line index and (5 — 1 . . .Q labels the six modes as- 
sociated with each cutting Une. The phonon wavevector 
q is defined on a one-dimensional Brillouin zone given by 
\q\ < vr/T. 

The second quantized phonon Hamiltonian is given by 

= ^/3^(9) bl^q bf3^q (3) 

where bp^^ creates a phonon in a state with energy 



transition to quantum mechanics and define a symmetry- 
adapted second-quantized phonon displacement operator 

UrJ = 5I^'3..(g) S{9j) er{l3,q,v) 

with the quantized phonon amplitudes 

2. ^MN-^^^ hujp^{q) 

Using the phonon displacement operator in Eq. ([6|) we 
arrive at the desired second-quantized electron-phonon 
Hamiltonian 

^-^p - EEE mA^^i) 

ss'P fJ.l' kq 

^ '^l' ,fj. + I^,k + q ^Sfj,k (^,3l/g + ^/3iy,-g) (8) 

The first term in the electron-phonon Hamilto- 
nian, Eq. ((HI), is proportional to the operator 
4',f,+i',k+q bfjuq- This term describes a phonon ab- 
sorption process in which an electron with energy Eg^ (k) 
absorbs a phonon with energy Huji3^{q) and then scat- 
ters into the state with energy Esi^^+^{k -\- q). Likewise, 
the second term in Eq. ([8]) describes a phonon emission 
process. 

Details concerning the evaluation of the electron- 
phonon interaction matrix elements Mj^'^, p{k,q) are 
given in Appendix [C] In evaluating these matrix ele- 
ments, we make explicit use of 2pz atomic wavefunctions 
and screened atomic potentials obtained from an ab initio 
calculation in graphene. 



C. Electron-phonon coupling 

The position-dependent single-electron Hamiltonian 
for a single-walled carbon nanotube is given by 

i?(r)=ro + ^Wc(r-(R.j+U,j)) (4) 

rJ 

where the kinetic energy is Tq and Wc(r) is the carbon 
atom potential. In Eq. ^ R^j is the equilibrium position 
of the r-th atom in the J-th two-atom unit cell and j is 
the corresponding atomic displacement from equilibrium 
as defined in Appendix [Bl 

Expanding this Hamiltonian in a Taylor series to first 
order in the atomic displacements, we obtain a position- 
dependent electron-phonon interaction Hamiltonian 

Hep{r) = - J2^v,{r- R,j) • U,j (5) 

rJ 

Starting from the classical symmetry-adapted atomic 
displacement in Eq. (|B7|) of Appendix [Bl we make the 



D. Coherent phonon generation 

The second quantized phonon displacement operators 
in Eq. ([6]) are defined in terms of a sum over phonon 
modes of the second quantized operator bj^j^^+bji^ For 
each phonon mode in the nanotube we are thus motivated 
to define a coherent phonon amplitude given hy^ 

Qpuq{t) = {bfi,q{t) +bl^_^{t)) 

where ( ) denotes the statistical average. Equations of 
motion for Qpuqit) can be obtained from the phonon and 
electron-phonon Hamiltonians in Eqs. ^ and ([5]). From 
the Heisenberg equation we obtain 

0^2 + ^/3'>(9) Qp'^qit) = 

ss' 

In coherent phonon spectroscopy, we assume that the 
optical pulse and hence the distribution of photoexcited 
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carriers is spatially uniform over the nanotube. In this 
case the electronic density matrix is diagonal and can be 
expressed as 

(cI',^+^,fe-gW c^pfcW) = 6,^,, Sq^o fsf,ik,t) (10) 

where fs^{k,t) is the electron distribution function in 
subband s/i with wavevector 

The only coherent phonon modes that are excited are 
the V — q — modes whose amplitudes satisfy a driven 
oscillator equation 



at2 



■iulQpit)^Spit) 



(11) 



where Qf}{t) = Qpoo{t) and ujp = Lop,^{(i = 0). There 
is no damping term in Eq. (jlip since anharmonic terms 
in the electron-phonon Hamiltonian are neglected. We 
solve the driven oscillator equation subject to the initial 
conditions (3/3(0) = 0/3(0) = 0. Taking the initial con- 
dition into accomit, the driving hmction Sfj{t) is given 

by 



sf^ik 



Xk,t)-f%{k)) (12) 



where fsfj.{k, t) are the time-dependent electron distribu- 
tion functions, f^s^Xk) are the initial equilibrium electron 
distribution functions, and M^^{k) = M^J^(fc, q — 0). 

The coherent phonon driving function Si3{t) depends 
on the photoexcited electron distribution functions. In 
principle, we could solve for the time-dependent distribu- 
tion functions in the Boltzmann equation formalism tak- 
ing photogeneration and relaxation effects into account. 
In CP spectroscopy, an ultrafast laser pulse generates 
electron-hole pairs on a time scale short in comparison 
with the coherent phonon period. In our experimental 
work, we typically use 50 fs ultrafast laser pulses to ex- 
cite RBM coherent phonons with oscillation periods of 
around 0.14 ps (fe^/j « 30 meV or 242 cm~^). 

After photoexcitation the electron-hole pairs slowly 
scatter and recombine. Jiang et aZ.— recently carried out 
a study of the clcctron-phonon interaction and relaxation 
time in graphite and found relaxation times on the order 
of a few picoseconds which is much slower than either 
the ultrafast laser pulse or the RBM coherent phonon 
oscillation period. The driving function Si3{t) thus rises 
sharply in a step-like fashion and then slowly vanishes as 
the distribution functions /s^(fc,t) return to fsp,ik)- The 
rapid initial jump in Sis{t) gives rise to an oscillatory part 
of the coherent phonon amplitude Qfjlt) at the coherent 
phonon frequency lo^ while the slow subsequent decay of 
Sp{t) gives rise to a slowly varying background. Since the 
observed CP signal is proportional to the power spectrum 
of the oscillatory part oiQp{t), we choose to ignore relax- 
ation effects and retain only the rapidly varying photo- 
generation term in the Boltzmann equation. Neglecting 
carrier relaxation will have a negligible effect on the com- 
puted CP signal since the relaxation time is much greater 
than the coherent phonon period. 



The photogeneration rate in the Boltzmann equation 
depends on the polarization of the incident ultrafast laser 
pulse. Using an effective mass model, Ajiki and Ando^ 
showed that optical absorption in an isolated single- 
walled carbon nanotube for polarization perpendicular 
to the nanotube axis is almost perfectly suppressed by 
photo-induced charge (the depolarization effect). Re- 
cently, Popov and Henrard*'' undertook a comparative 
study of the optical properties of carbon nanotubes in 
orthogonal and nonorthogonal tight-binding models and 
found that optical absorption due to light polarized par- 
allel to the tube axis (z axis) is greater than absorption 
due to light polarized perpendicular to the axis by about 
a factor of five. Consequently we confine our attention to 
light polarized parallel to the tube axis only. We compute 
the photogeneration rate in the electric dipole approxi- 
mation using Fermi's golden rule. In the case of parallel 
polarization, optical transitions can only occur between 
states with the same angular momentum quantum num- 
ber jj.J^'^ For the photogeneration rate we find 



dt 



uit) 



gen 

■ [fs'^.{k,t)- 



!,^{k,t)) s{^E':,,{k)-hJ) (13) 



where AiJ^'^,, (fc) = |i?s^(fc) — £'s'p(fc)l are the k dependent 
transition energies, hw is the pump energy, u{t) is the 
time-dependent energy density of the pump pulse, e is 
the electron charge, mo is the free electron mass, and 
Ug is the index of refraction in the surrounding medium. 
The optical matrix element is given by 

p!:Ak)^^Y.c*As',t,,k) 

v2mo ^ — ^ 

V ^ j.t 

C^{s,^i,k) e*"^^('=^'') M^(r',rJ) (14) 

rJ 

where the sum over r J is taken over fourth nearest neigh- 
bors of the atom at Rr'o. In Eq. (fT4|) . C^(s, /i, k) are the 
expansion coefficients for the symmetry- adapted ETB 
wavefunctions obtained by solving the matrix eigenvalue 
problem in Eq. (|A8p of Appendix YK\ and (j)j{k,ijL) is the 
phase factor defined in Eq. (jA7[) . 

The z components of the atomic dipole matrix ele- 
ments (which can be evaluated analytically) are given by 



f d 

^ dr ip*,o{r - R^'o) ^ <Prj(r - R-rj) 



(15) 



where the 2pz orbitals (prj are defined in Eq. (jC3|) . Note 
that the squared optical matrix element \Pss'ik)\'^ has 
units of energy. We point out that optical dipole matrix 
elements in the vicinity of the K point in both graphite 
and carbon nanotubes have been studied previously in 
Ref. SH. 

The pump energy density u{t) is related to the fluence 
F — J dt u{t) {c/ng). To simplify our theoretical model, 
it is assumed that the pump beam consists of a train 
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of -^puise identical Gaussian pulses each with an inten- 
sity full width at half maximum (FWHM) of Tp and a 
Lorentzian spectral lineshape with a FWHM of Tp. The 
Gaussian pulses are equally spaced in time with the time 
interval between pulses being T^p^igg- The peak intensity 
of the first pulse is taken to occur at t = 0. To account 
for spectral broadening of the laser pulses we replace the 
delta function in Eq. (fT^ with'*^ 



S{AE - hw) 



Tp/{27r) 



(AE - hiof + (rp/2)2 



(16) 



From the coherent phonon amplitudes, the time- 
dependent macroscopic displacements of each carbon 
atom in the nanotube can be obtained by averaging 
Eq. dH). Thus 

U.j(t) =Y.^P ^(^J) Qpit) (17) 


where = A;3,o(0) and ef = er(/3, 0,0). 

It is apparent that only four coherent phonon modes 
can be excited in a carbon nanotube regardless of the 
chirality. Since v = q = the CP active mode frequen- 
cies and polarization vectors are found by diagonalizing 
a single 6x6 dynamical matrix in Eq. (jB10|) . Two of the 
six mode frequencies top are zero and these phonons are 
not excited since the driving term Si3{t) vanishes as can 
be seen in Eq. (H^). Of the remaining four modes, the 
one with the lowest energy is the radial breathing mode. 

Coherent acoustic phonon modes whose energies van- 
ish at q — cannot be excited in an infinitely long car- 
bon nanotube under conditions of uniform illumination 
by the pulse laser. If, however, the electric field of the 
pump laser could be made to vary spatially along the 
nanotube axis with a periodicity given by a real space 
wavevector ^Zp^^gg, it would be possible to generate co- 
herent acoustic phonons which would travel along the 
nanotube at the acoustic sound speed. The generation 
of coherent acoustic phonons has been demonstrated in 
semiconductor superlattices where the pump laser gen- 
erates carriers in the quantum wells thus giving rise to 
carrier distribution functions having the periodicity of 
the superlattice4^i^i^i^ In the case of coherent acoustic 
phonons in superlattices, the coherent phonon lattice dis- 
placement satisfies a driven loaded string equation rather 
than a driven harmonic oscillator equation.'^ 



E. Absorption spectrum 

In coherent phonon spectroscopy a probe pulse is used 
to measure the time- varying absorption coefficient of the 
carbon nanotube. The time-dependent absorption coef- 
ficient is given by^i^ 



where e2{fu^,t) is the imaginary part of the time- 
dependent dielectric function evaluated at the probe pho- 
ton energy huj. 

The imaginary part of the nanotube dielectric function 
is obtained from Fermi's golden rule 



STT^e^ f h^\^ f dk 



E 



a{huj, t) 



huj 

UghC 



(18) 



(/.mW-A'mW) 5(^E,,^{k)~E,^{k)-hio^ (19) 

where At = Tr{dt/2)'^ is the cross-sectional area of the 
tube and dt is the nanotube diameter. In our model, 
we replace the delta function in Eq. (|19p with a broad- 
ened Lorentzian spectral lineshape with a FWHM of 
Ts- The distribution function fsii{k) and bandstruc- 
ture Es^(k) are time-dependent. The time-dependence 
of fsfj,{k) comes from the Boltzmann carrier dynamics 
which can include the photogeneration of electron-hole 
pairs as well as various carrier relaxation effects. The 
time-dependence of Es^{k) arises from variations in the 
carbon-carbon bond lengths due to the coherent phonon 
induced atomic displacements Urj(t) given in Eq. (jl7|l . 
This time-dependent deformation of the nanotube bond 
lengths alters the tight-binding Hamiltonian and overlap 
matrix elements in the extended tight-binding model de- 
scribed in Section IIII Al Note that to first order in the 
lattice displacements the energies Es^ (fc) vary with time 
while the tight-binding wavefunctions and optical matrix 
elements Pgg,{k) do not. 



F. Coherent phonon spectrum 

In the coherent phonon spectroscopy experiments de- 
scribed by Lim et al. in Refs. [H and [ij single color 
pump-probe experiments are performed on an ensemble 
of nanotubes. The excitation of coherent phonons by the 
pump modulates the optical properties of the nanotubes 
and gives rise to a transient differential transmission sig- 
nal. After subtraction of a slowly varying background 
component, the coherent phonon spectrum is obtained 
by taking the power spectrum of the time-dependent dif- 
ferential transmission. In our model we simulate single- 
color pump-probe experiments and take the theoretical 
CP signal to be proportional to the power spectrum of 
the transient differential transmission after background 
subtraction. We compute the power spectrum using the 
Lomb periodogram algorithm described in Ref. IZ^. We 
find that using the Lomb periodogram to evaluate the 
power spectrum in our theoretical results is more con- 
venient than using fast fourier transform methods since 
it works well for data sets whose size is not an integer 
power of two or whose data points are not evenly spaced. 

It is worth emphasizing that in CP spectroscopy 
we measure the power spectrum of the time-dependent 
coherent phonon-modulated differential transmission. 
Thus, as a function of the pump energy, the CP spec- 
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k (ti/T) 

FIG. 8: (color online) The upper panel show the four lowest 
lying electronic tt bands in the one-dimensional Brillouin zone 
for an (11,0) zigzag nanotube in the extended tight-binding 
(ETB) model. All the bands are doubly degenerate. The 
allowed optical transition Ea for light polarized along the 
tube axis are indicated. The lower panel shows the square of 
the optical matrix elements defined in Eq. [14] for each of the 
allowed optical transition. 



trum at the coherent phonon frequency tracks the abso- 
lute value of the first derivative of the static absorption 
coefficient. This is nicely illustrated in Fig. 6 of Ref. [l^ 
where it is shown that an excitonic peak in the absorp- 
tion spectrum will give rise to a symmetric double peaked 
structure in the CP power spectrum. 



IV. THEORETICAL RESULTS 

To illustrate our theoretical model, we will discuss in 
some detail simulated CP spectroscopy experiments in 
an undoped (11,0) zigzag nanotube. This is a mod 2 
semiconducting nanotube with mod (n — m, 3) = 2 be- 
longing to the family of nanotubes with 2n + m — 22. We 
choose this example since we have performed CP spec- 
troscopy for the (11,0) nanotube (see Fig. [TH]) and be- 
cause Lim et al. measured coherent lattice vibrations in 
a micelle-suspended solution of carbon nanotubes with a 
diameter range of 0.7 to 1.3 n m^^i^^ and found strong CP 
signals due to excitation of coherent REM modes in this 



family of nanotubes. 



A. Bandstructure and absorption spectra 

The four lowest-lying one-dimensional electronic tt 
bands for the (11,0) nanotube are shown in the upper 
panel of Fig. [51 In our model we ignore structural op- 
timization and assume that the carbon atoms lie on the 
surface of the rolled up unreconstructed graphene cylin- 
der. In an unreconstructed zigzag nanotube the length of 
the translational unit cell is T = Vs a where a = 2.49 A 
is the hexagonal lattice constant in graphene.- In Fig. [5| 
the conduction bands have positive energy and the va- 
lence bands have negative energy. Since the electronic 
problem has been reduced to solving the 2x2 matrix 
eigenvalue problem in Eq. (|A8|) . the conduction bands in 
Fig. [8] with a given value of the angular quantum number 
/i can only mix with the valence band having the same 
value of /X. The four bands shown in the figure are dou- 
bly degenerate with two distinct values of /i giving rise to 
the same band energies. The allowed optical transitions 
for z-polarized light (with selection rule A/u = 0) are 
indicated by vertical arrows and are labeled En . . . -B44. 

The lower panel of Fig. [5| shows the square of the op- 
tical matrix elements defined in Eq. ([M]) for each of the 
transitions En.. .£^44. For the -En, E22 and -B33 tran- 
sitions the squared optical matrix elements are strongly 
peaked at the van Hove singularity at the direct band 
gaps. The size of this peak in the squared optical matrix 
elements increases as the band gap moves away from the 
Fermi energy at E = 0. The absorption for these three 
transitions are sharply peaked at the band edge due to 
the van Hove singularity in the joint density of states as 
well as the peak in the squared optical matrix elements 
that occurs there. For the ,£44 transition, the conduction 
and valence bands are very flat giving rise to an enhanced 
van Hove singularity in the joint density of states while 
the squared optical matrix element is a slowly varying 
function of k. For this transition, the peak in the ab- 
sorption spectrum is due almost entirely to the sharply 
peaked joint density of states. 

As a check on our theoretically calculated squared opti- 
cal matrix elements, we compared our results with optical 
dipole matrix elements calculated independently by Jiang 
et al. in Ref. [4^| for the metallic (5,5) armchair and (6,0) 
zigzag tubes for light polarized parallel to the tube axis. 
Our squared optical matrix elements are proportional to 
the square of the dipole matrix elements shown in Figs. 2 
and 3 in Ref. |49|. We found excellent agreement between 
our theoretical results and the corresponding results for 
the two tubes considered. 

With the electronic band structure and squared op- 
tical matrix elements shown in Fig. [51 we can obtain 
the absorption coefRcient of the (11,0) nanotube using 
Eqs. P5)l and and the carrier distribution functions. 
The computed absorption coefficient for the undoped 
(11,0) nanotube in thermal equilibrium at room tem- 
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FIG. 9: Computed absorption coefficient as a function of 
photon energy for all nanotubes (n, m) in the family for which 
2n + m = 22. The photons are assumed to be linearly polar- 
ized with the electric polarization vector parallel to the tube 
a:xis. Exciton effects are neglected. 
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FIG. 10; (color online) Computed photoexcited carrier den- 
sity per unit length for nanotubes in Family 22 after pho- 
toexcitation by a 50 fs Gaussian laser pulse with a fluence of 
10"^ J/cm^ pumping at the peak of the E22 features in Fig.|9] 



perature is shown in the upper curve of Fig. [9] where the 
spectral FWHM Hnewidth is taken to be Ts = 0.15 eV. 
Also shown in the figure are the absorption spectra of the 
other members of Family 22, namely the (10,2), (9,4) 
and (8, 6) nanotubes. 



B. Generation of coherent RBM phonons 

In a typical simulation, we excite coherent RBM 
phonons with a single 50 fs Gaussian laser pulse pump- 
ing at the peak of the broadened E22 transitions shown in 
Fig. m The pump fluence is taken to be 10^^ J/cm^, the 
FWHM spectral linewidth is assumed to be Fp = 0.15 eV, 
and the time scale is chosen so that the pump reaches its 
peak intensity at t = 0. The pump energies huj at the £^22 
peaks are taken to be 2.05 eV for the (11, 0) nanotube, 
2.04 eV for (10,2), 1.97 eV for (9,4), and 1.89 eV for 
(8, 6). The photogenerated carrier distribution functions 
are obtained from Eq. (|13p and the photoexcited carrier 
densities per unit length are shown in Fig. [10] for Family 
22 nanotubes. The carrier densities after photoexcitation 
all lie in the range from 70 to 90 cm~^ increasing as we 
go from (11,0) to (8,6). 

Using the photogenerated carrier densities fsfj.{k, t) we 
can obtain the coherent phonon amplitudes Qfj{t) by 
solving the equation of motion (jlip with the driving func- 
tion Sp{t) given in Eq. (fT^. As we noted earlier, (3 la- 
bels the six coherent phonon modes for each value of /i. 
Coherent phonon oscillations can only be excited in the 
four q — Q modes with non-zero frequency corresponding 
to /X = 0. The six /i = phonon dispersion curves are 




FIG. 11: Phonon dispersion relations for fi — phonons in 
(11,0) nanotubes. Only the four q — modes with non-zero 
frequency can be excited in coherent phonon spectroscopy. 



shown in Fig.[TT]for the (11, 0) nanotube. At g = 0, there 
are two acoustic modes with zero frequency. The one 
with the lower sound speed is the twisting mode (TW) 
in which the A and B atom sublattices move in phase in 
the circumferential direction. The mode with the higher 
sound speed is the longitudinal acoustic mode (LA) in 
which the A and B atom sublattices move in phase along 
the tube axis. The remaining four q = modes are the 
CP active modes. The lowest CP active mode is the ra- 
dial breathing mode at 37.1 meV (300 cm"^) in which all 
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FIG. 12: (color online) Driving function kernel for RBM 
coherent phonons in the (11,0) nanotube for several optical 
transitions En in the impulsive excitation approximation. 



the atoms move in phase radially in and out. The next 
highest CP active mode is the out of plane transverse 
optical mode (oTO) at 117 meV (944 cm~^) in which 
the A atom sublattice moves radially outward as the B 
atom sublattice moves radially inward. There are two 
closely spaced high frequency modes. The lower of these 
at 194.8 meV (1571 cm~^) is the in-plane transverse op- 
tical mode (iTO) in which the A and B atom sublattices 
move in the circumferential direction out of phase with 
each other. The highest mode is the longitudinal opti- 
cal mode (LO) at 197.9 meV (1596 cm~^) in which the 
A and B atom sublattices move out of phase with each 
other along the tube axis. 

In our model we neglect slow carrier relaxation ef- 
fects and retain only the photogeneration term in the 
Boltzmann equation. The net photogenerated conduc- 
tion band electron distribution function /cp(fc) — /c^(fc) 
is then equal to the net photogenerated hole distribution 
function for each value oi k. In this case, we can obtain 
a simplified expression for the coherent phonon driving 
function which only involves the photogenerated conduc- 
tion band electron distributions. We find 

Spit) = ^^(^) (/^m(^) - (20) 

fj.k 

The driving function kernel S^{k) is given by 

= (Mf^(fc) - M^P^ik)) (21) 

where Mf (fc) (s = v,c) are the same matrix elements 



appearing in Eq. ([12]). Each value of fi in the impul- 
sive excitation model corresponds to a specific optical 
transition En. The k dependence of Sf^{k) for the RBM 
phonon in the (11,0) nanotube is shown in Fig. [T^ for 
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FIG. 13: (color online) Coherent phonon generation in (11,0) 
nanotubes by photoexcitation at the En and E22 transition 
energies. The upper panel shows the density of photoexcited 
electron-hole pairs per unit length, the middle panel shows 
the coherent phonon driving function, and the bottom panel 
shows the RBM coherent phonon amplitude. 



the first four optical transitions. For the RBM mode, we 
chose the unit mode polarization vector to point radially 
outward so that positive values of S{k) contribute to a 
radially outward directed driving term. 

As can be seen in Fig. [T^] both positive and negative 
values of S{k) are possible. If, for example, we were to 
pump near the En band edge, the electron distribution 
functions would be localized near fc = and we would get 
negative values for S{t). The signs of the driving function 
kernels near fc = for the En transitions are in agreement 
with other results reported in the literature?^ The sign 
of (k) is the negative of the sign of the electron-phonon 
matrix element M^^{k) — M^^{k) appearing in Eq. (|2ip . 
This matrix element has been obtained by Machon et 
al. in an ab initio calculation reported in Ref. for 
the lowest four optical transitions with light polarized 
parallel to the tube axis. In Table I of Machon et ali^ 
the sign of the band edge electron-phonon matrix element 
for Ell in the (11,0) tube has a sign opposite to that of 
the higher lying transitions in agreement with the results 
shown in Fig[T2l 

In Fig.[T3]we plot the photoexcited carrier density n{t), 
the coherent phonon driving function S{t), and the coher- 
ent phonon amplitude Q{t) for RBM coherent phonons 
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the E22 transition energy. The dimensionless coherent 
phonon ampUtudes Q{t) for the radial breathing mode 
obtained by solving the driven oscillator equation (jlip 
are shown in the lower panel of Fig. [T3] where the curves 
are offset for clarity. As expected, the coherent phonon 
amplitudes oscillate at the RBM frequencies about a new 
positive equilibrium point. We note that for the radial 
breathing mode the coherent phonon amplitude is pro- 
portional to the differential change in the tube diameter. 

In Fig. [141 we see that the magnitude of the coher- 
ent phonon oscillations depends on chirality. The size 
of the jump in the driving function Sit) due to pho- 
toexcitation, and therefore the magnitude of the oscil- 
lations in Q{t), should be roughly proportional to the 
product of the electron-photon and electron-phonon ma- 
trix elements. The electron-phonon interaction matrix 
element for E22 transitions in mod 2 tubes is largest 
in the zigzag nanotube limit while the corresponding 
electron-photon interaction matrix element is largest in 
the armchair limit 4^ Thus we expect S{t) and Q{t) to 
have maxima somewhere between these two limits. 



C. Resonant excitation of coherent phonons 



FIG. 14: (color online) Upper panel: Coherent phonon driv- 
ing functions for RBM oscillations in Family 22 nanotubes 
photoexcited by 50 fs pumps at the broadened E22 absorp- 
tion peaks in Fig. [O] Lower panel: Corresponding coherent 
phonon amplitudes. Curves are offset for clarity. 



in (11,0) tubes for 50 fs z-polarized laser pulses with 
photoexcitation energies of 1.07 eV and 2.05 eV. These 
correspond to the En and E22 absorption peaks seen in 
Fig. [5] respectively. The absorption peaks are compara- 
ble for £'11 and i?22- However, for constant fluence, the 
number of photoexcited carriers n cx a/hw (see Ref.l42, p. 
341) and so the number of photoexcited carriers is larger 
for the Ell transition primarily as a result of the smaller 
transition energy. As expected from Fig. [121 the coher- 
ent phonon driving functions Sit) and amplitudes Q{t) 
have different signs in the two cases. This means that 
for photoexcitation at the En transition in mod 2 nan- 
otubes the tube diameter decreases and oscillates about 
a smaller equilibrium diameter while the opposite is true 
for photoexcitation at the £^22 transition energy. As an 
aside, we note that in mod 1 tubes the predicted E22 
and £11 photoexcited diameter oscillations have the op- 
posite initial phase relative to the corresponding mod 2 
oscillations as will be show in Section IIV El 

If we pump at the broadened E22 absorption peaks 
with 50 fs pulses of the same spectral width for each 
nanotube in Family 22, we obtain the driving functions 
shown in the upper panel of Fig. [141 In all cases, the 
driving functions are positive which implies that all the 
nanotubes in Family 22 initially move radially outward 
in response to photoexcitation by an ultrafast pump at 



In our experimental work, we resonantly excite coher- 
ent RBM phonons in specific chirality nanotubes in a mi- 
celle suspended sample by using a train of pump pulses 
with a repetition rate equal to the RBM period. To il- 
lustrate resonant excitation of coherent phonons in nan- 
otubes, we repeated our simulations for the (11,0) nan- 
otube using a train of six Gaussian pulses pumping at 
the E22 transition energy (2.05 eV) where the pulse rep- 
etition period was in phase and 180 degrees out of phase 
with the period of the RBM phonon. 

The results are illustrated in Fig. 1151 The upper panel 
of Fig. [T51 shows the pump laser intensity as a function of 
time for the in-phase and out-of-phase cases. For compar- 
ison, we also plot the pump laser intensity for the single 
50 fs Gaussian pulse used in our earlier simulation. In 
all three cases the fluence is taken to be 10~^J/cm^ and 
hence the final density of photogenerated carriers were 
the same. We note that in the figure, the intensities of 
the two pulse trains are multiplied by a factor of six since 
the intensities of the Gaussian pulses scale inversely with 
the number of pulses if the fluence is held constant. 

The corresponding coherent phonon amplitudes for the 
RBM phonons are shown in the lower panel of Fig. [THl 
For the single Gaussian pulse, the coherent phonon am- 
plitude is the same as that shown in Fig.[T31for the (11,0) 
nanotube. For the in-phase case, the coherent phonon 
amplitude is magnified after each Gaussian pulse and at 
long times is identical to the coherent phonon amplitude 
obtained using the single Gaussian pulse. Since it is the 
long time behavior of the oscillating part of the coherent 
phonon amplitude that determines the GP signal, the 
two resulting GP spectra are also identical. When the 
Gaussian pulses are 180 degrees out of phase, the oscil- 
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FIG. 15: (color online) Upper panel: Pump laser intensity vs. 
time for a single 50 fs Gaussian pulse and two trains of six 50 
fs Gaussian pulses with repetition periods in phase and 180 
degrees out of phase with the (11,0) RBM phonon period. 
In all cases, we pump at the (11,0) E22 transition energy 
(2.05 eV) with the same fluence. Note that the intensities of 
the two pulse trains are magnified by a factor of six relative 
to the intensity of a single Guassian pulse. Lower panel: Cor- 
responding coherent phonon amplitudes for the (11,0) RBM 
coherent phonon oscillations. 



lating part of the coherent phonon amphtude and hence 
the resulting CP spectra are completely suppressed. 

In Fig. [15] the predicted coherent phonon amphtudes at 
long times oscillate rapidly about the steady state value 
Qf3 ~ Sf){t 00) /ujp. In real nanotubes the driving 
function Sf3{t) will slowly vanish as the carriers recom- 
bine. In addition the coherent phonon amplitudes will 
slowly decay on a time scale of tens of picoseconds as 
evidenced in Fig. [31 



D. Coherent phonon detection in mod 2 tubes 

The generation of coherent phonons results in periodic 
oscillations of the carbon atoms which in turn modulate 
the optical properties of the nanotube. These coherent 
phonon oscillations can be detected by measuring tran- 
sient optical properties in pump-probe experiments. We 
simulate single color pump-probe measurements in which 
we pump with light linearly polarized along the tube axis 
and measure the transient differential gain as a function 
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FIG. 16: (color online) The lower panel shows the absorption 
coefficient in (11,0) nanotubes as a function of photon energy 
for light polarized along the tube axis. The upper panel is 
a contour map showing the coherent phonon power spectrum 
as a function of pump-probe energy and phonon frequency. 



of probe delay for a probe pulse having the same en- 
ergy and polarization as the pump. The coherent phonon 
spectrum is obtained by scanning the pump-probe energy 
and in our simulations, the pump fluence, duration, and 
FWHM spectral linewidth are assumed to be constant as 
we vary the pump-probe energy. 

The coherent phonon (CP) spectrum for mod 2 (11, 0) 
nanotubes is shown in Fig. 1161 The bottom panel shows 
the absorption spectrum for light linearly polarized along 
the tube axis assuming a FWHM linewidth of 0.15 eV. 
The upper panel shows a contour map of the coherent 
phonon power spectrum as a function of pump-probe en- 
ergy and photon energy. The CP intensity is proportional 
to the power spectrum and as we scan in photon energy 
two large peaks are observed at the En and E22 tran- 
sitions at a phonon energy near 37.1 meV (300 cm~^) 
which corresponds to the RBM frequency of the (11,0) 
nanotube. Comparing the upper and lower panels of 
Fig. 1161 we can verify our earlier assertion that as we 
scan in pump energy the CP intensity at the RBM co- 
herent phonon frequency is proportional to the absolute 
value of the first derivative of the absorption coefficient. 
We should point out that our theoretical model does not 
include many-body Coulomb effects and so the position 
and shape of the CP signal is due to modulation of the 
free carrier En and E22 transitions. Thus in our free 
carrier model, we see an asymmetric double peak at each 
transition with the stronger peak at low pump energy 
and the weaker peak at higher energy. 

In our simulations with 50 fs laser pulses, we excite 
coherent RBM phonon modes since the pulse duration is 
much less than the RBM phonon oscillation period. To 
excite the higher lying coherent phonon modes it is nec- 
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FIG. 18: (color online) Coherent phonon intensity at the 
RBM frequency as a function of pump-probe energy for sev- 
eral mod 2 semiconducting nanotubes at the E22 transition. 
The experimental CP spectra are in the right panel and the 
simulated CP spectra are in the left panel. The upper four 
curves in each panel are for nanotubes in Family 22 and the 
lower four curves are for tubes in Family 25. Each curve is 
labeled with the chirality (n, m) and the RBM phonon energy 
in meV and offset for clarity. 



FIG. 17: (color online) Coherent phonon spectra for coherent 
RBM and LO phonons in (11,0) nanotubes photoexcited by 
an ultrafast z-polarized 5 fs pump. The CP spectra at the 
RBM and LO energies (37.1 and 198 meV) are plotted as a 
function of pump-probe energy. Two features are seen near 
Ell and £^22 transitions. The LO curve is multiplied by 1000 
and offset for clarity. 



essary to use shorter laser pulses. Recently, Gambetta et 
al}^ used sub-lO-fs laser pulses to excite the RBM and 
LO modes vifhile excitation of the oTO and ITO modes 
were not observed. To examine this case, we simulated 
CP spectroscopy in (11,0) nanotubes using short 5 fs 
laser pulses. In qualitative agreement with the measure- 
ments of Gambetta et ai, we find that the two strongest 
modes are the RBM and LO modes while the strengths 
of the oTO and iTO modes are found to be negligible. 
Our result for the iTO mode is consistent with the chiral- 
ity dependent Raman G-band intensity in which the iTO 
signal is absent in zigzag nanotubes .^fifi The CP spectra 
for the RBM and LO modes are shown in Fig. [T71 where 
the differential gain power spectra at the RBM and LO 
energies (37.1 and 198 meV, respectively) are plotted as a 
function of pump-probe energy. Two strong features are 
seen near the En and E22 transition energies. In this 
example, the two curves have similar shapes but the LO 
CP spectrum (multiplied by a factor of 1000 in the fig- 
ure) is much weaker than the corresponding RBM spec- 
trum. We note that shortening the duration of the laser 



pulse from 50 fs to 5 fs enhances the strength of the E22 
peak relative to the En peak. This can best be seen by 
comparing the bottom curves in Fig. 1171 with the bottom 
curve in Fig. 1211 

It is useful to examine trends in the CP spectra within 
and between mod 2 semiconducting nanotube families by 
plotting the theoretical CP intensity at the RBM phonon 
frequency as a function of pump-probe energy. This is 
done in the left panel of Fig. [TH] where we plot our theo- 
retical CP intensity at the RBM frequency as a function 
of pump-probe energy for all nanotubes in Families 22 
and 25. The curves for each nanotube are labeled with 
the nanotube chirality (n, m) and the RBM phonon en- 
ergy in meV. In each nanotube, we see peaks in the CP 
spectra corresponding to E22 transitions. Within a given 
family, the CP intensity tends to decrease as the chi- 
ral angle increases, i.e., as the chirality goes from (n, 0) 
zigzag tubes to (n, n) armchair tubes. From Fig. [18] we 
can also see that the theoretical CP intensity increases 
as we go from Family 22 to Family 25. 

The right panel of Fig. [TS] shows the corresponding ex- 
perimental CP spectra for the nanotubes in Families 22 
and 25. Comparing experimental and theoretical curves 
in Fig. [THl we see that our theory correctly predicts the 
overall trends in the CP intensities. Since we are using 
pump probe methods to study an ensemble of micelle- 
suspended nanotubes, the relative agreement between 
the theoretically calculated and experimentally measured 
CP intensities suggests that nanotubes of different chi- 
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ralities in Families 22 and 25 in the micelle-suspended 
sample studied are equally probable and that the mea- 
sured CP signal strengths are an intrinsic property of the 
tubes. 

There are discrepancies in the predicted pump-probe 
energies of the peaks on the order of 0.4 eV or less. A 
discrepancy of this size is expected since we have not 
included many-body Coulomb interactions in our the- 
oretical model. It is well established that both the 
excitonic red shift and the self-energy blue shift are 
very large in nanotubes, with the latter exceeding the 
former i^i^iSiSi^i^ We also note that the dielectric func- 
tion of the surrounding medium also influences the exci- 
tonic transition energiesi^ There are also differences be- 
tween the theoretical and experimental CP lineshapes in 
Fig. [iHl If we compare theoretical and experimental CP 
spectra for the (12,1) nanotube, we see that both exhibit 
a double peaked structure. However, the lower energy 
theoretical peak is much stronger than the higher energy 
peak whereas the two experimental peaks have compa- 
rable strength. This discrepancy can be attributed to 
strong excitonic modification of the shape of the nan- 
otube absorption spectrum whose time-dependent mod- 
ulation of the probe pulse gives rise to the shape of the 
CP signal. The free carrier absorption edge is highly 
asymmetric while the excitonic absorption spectrum ex- 
hibits a symmetric peak at the exciton transition energy, 
thus accounting for the discrepancy. Our theory qualita- 
tivley agrees with experiment, but to obtain quantitative 
agreement, one must include details of the Coulomb in- 
teraction. 

Our experimental and theoretical results are also in 
qualitative agreement with the results of CP spectroscopy 
measurements previously reported by Lim et al. in 
Ref. [13. As these authors note, the tendency of the CP 
intensity to increase with family index is in contrast to 
the situation in resonant Raman scattering where the 
strength of the resonant Raman signal is observed to de- 
crease as the family index increases. 



E. Coherent phonon detection in mod 1 tubes 

It is useful to perform a comparison between mod 1 
and mod 2 semiconducting nanotubes with the same 
chiral angle. To this end, we compare the Family 22 
(11,0) mod 2 nanotubes with Family 26 (13,0) mod 1 
tubes. In both cases, the pump lasers have the same 
fluence (10~^ J/cm^), pulse duration (50 fs), and spec- 
tral linewidth (0.15 eV), while the pump energies cor- 
respond to the maxima in the broadened E22 absorp- 
tion features. For the (11,0) tube we pump at 2.05 eV 
and for the (13,0) tube we pump at 1.84 eV. The 
time-dependent photogenerated carrier densities per unit 
length are shown in the upper panel of Fig. 1191 The mid- 
dle and lower panels show the coherent phonon driving 
functions and corresponding coherent phonon amplitudes 
for REM coherent phonons in the two cases. For (11,0) 
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FIG. 19: (color online) Coherent phonon generation in (11, 0) 
mod 2 and (13, 0) mod 1 semiconducting zigzag nanotubes by 
photoexcitation at the E21 transition energy. The upper panel 
shows the density of photoexcited electron-hole pairs per unit 
length, the middle panel shows the coherent phonon driv- 
ing function, and the bottom panel shows the RBM coherent 
phonon amplitude. 



mod 2 tubes, the coherent phonon driving function is 'pos- 
itive while for (13,0) mod 1 tubes, the driving function 
is found to be negative. Some insight into this behavior 
can be obtained by examining the /c-dependent driving 
function kernel for coherent RBM phonons in (13, 0) nan- 
otubes shown in Fig. [201 Comparing Fig. [201 with Fig.fT^ 
we see that the driving function kernels for En and E22 
transitions near fc = have opposite signs which accounts 
for the sign change. 

These results are supported by other studies reported 
in the literature^i^ The difference in the sign of the E22 
driving function kernels at fc = for the zigzag mod 2 
(11,0) tube shown in Fig. [T^ and the zigzag mod 1 (13,0) 
tube shown in Fig. [201 is due to a change in sign of the 
electron-phonon matrix element. This sign change in the 
electron-phonon matrix element for tubes with different 
mod numbers was also obtained independently by Jiang 
et al. in Ref. [H. Figure 1(a) in Jiang et ali^ shows 
the electron-phonon matrix element for coherent RBM 
phonons excited at the E22 transition as a function of 
the chiral angle 9 for mod 1 and mod 2 tubes (the SII 
and SI curves in Jiang's Fig. 1(a)). For zigzag tubes. 
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FIG. 20: (color online) Driving function kernel for RBM 
coherent phonons in the (13,0) nanotube for several optical 
transitions En in the impulsive excitation approximation. 



~ Q and from Jiang's Fig. 1(a) we can see that the 
sign of the electron-phonon matrix element for the mod 2 
(SI) tube is positive while the sign of the electron-phonon 
matrix element for the mod 1 (SII) tube is negative. This 
is consistent with the theoretical results shown in our 
Figs. [12] and [201 

We find that the CP intensity is very sensitive to the 
nanotube mod number. In general, the CP intensity in 
mod 2 semiconducting nanotubes is much larger than 
the CP intensity in mod 1 semiconducting tubes. This is 
illustrated in Fig. [21] where we plot the CP power spec- 
tra as a function of pump-probe energy at the RBM fre- 
quencies for the zigzag (13,0) mod 1 and (11,0) mod 2 
semiconducting nanotubes. For the (13,0) and (11,0) 
tubes the RBM phonon energies are 31.5 (254 cm~^) and 
37.1 meV (300 cm'^), respectively. Note that in Fig. [211 
the mod 1 curve is multiplied by a factor of 5. In general, 
we find that CP intensities in mod 2 tubes are consider- 
ably larger than CP intensities in mod 1 tubes. This is 
consistent with the experimental results of Lim et al. as 
reported in Refs.[l2 andfisl. We also find that in mod 1 
tubes, the En feature is more pronounced than the £'22 
feature in contrast to what is seen in the mod 2 case. 

The CP intensities as a function of chirality for nan- 
otubes in two mod 1 families (Families 26 and 29) are 
shown in Fig. [22] We find that in all cases the En fea- 
tures in the mod 1 tubes are much stronger than the £^22 
features. Within a mod 1 family, the CP intensity of the 
Ell feature is found to decrease with increasing chiral 
angle. We also find that the CP intensities decrease as 
the mod 1 family index increases. 
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FIG. 21: (color online) Coherent phonon power spectra at the 
RBM energies in the (13, 0) mod 1 and (11, 0) mod 2 semicon- 
ducting nanotubes as a function of pump-probe energy. The 
mod 1 CP power spectrum is multiplied by a factor of 5 and 
offset for clarity. 



V. SUMMARY 

Using femtosecond pump-probe spectroscopy with 
pulse shaping techniques, we have generated and de- 
tected coherent phonons in chirality-specific semicon- 
ducting single-walled carbon nanotubes. The signals 
were resonantly enhanced when the pump photon energy 
coincides with an interband exciton resonance, and anal- 
ysis of such data provided a wealth of information on the 
chirality-dependence of light absorption, phonon genera- 
tion, and phonon-induced band structure modulations. 

To explain our experimental results qualitatively and 
quantitatively, we have developed a microscopic theory 
for the generation and detection of coherent phonons in 
semiconducting single-walled carbon nanotubes via co- 
herent phonon spectroscopy. For extremely short laser 
pulses, we find that the two strongest coherent phonon 
modes are the RBM and LO modes. The CP spectrum 
of the LO mode is similar in shape to that of the RBM 
mode but is found to be much weaker. 

For the RBM modes, the CP intensity within a mod 2 
family tends to decrease with chiral angle, and the de- 
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FIG. 22: (color online) Coherent phonon intensity at the 
RBM frequency as a function of pump-probe energy for sev- 
eral mod 1 semiconducting nanotubes. The upper five curves 
are for nanotubes in Family 26 and the lower five curves are 
for tubes in Family 29. Each curve is labeled with the chiral- 
ity (n, m) and the RBM phonon energy in meV. The curves 
are offset for clarity. 



crease in CP intensity with chiral angle is found to be 
much more pronounced for the En feature. We also find 
that CP intensities are considerably weaker in mod 1 fam- 
ilies in comparison with mod 2 families. In general, the 
E22 CP intensities in mod 2 families are stronger than 
the Ell features. For RBM modes in mod 1 tubes, the 
£"11 intensities are stronger than the E22 intensities and 
tend to decrease with increasing chiral angle within a 
given family. As the family index increases, the En CP 
intensity in mod 1 tubes decreases. 

For mod 2 nanotubes, we predict that the tube diam- 
eter will initially increase for E22 photoexcitation and 
decrease for £^11 photoexcitation. In mod 1 nanotubes, 
the opposite is precdicted to be the case, i.e. the tube 
diameter will initially decrease for £^22 photoexcitation 
and increase for En photoexcitation. 

We compare our theoretical results with experimental 
CP spectra in mod 2 nanotubes and find that our theo- 
retical model correctly predicts the overall trends in the 
relative strengths of the CP signal both within and be- 
tween mod 2 families. We find discrepancies between our 



theoretical predictions with regard to the peak positions 
and lineshapes, which we attribute to Coulomb interac- 
tion effects that are not included in our calculations. 

For (11,3) mod 2 nanotubes, we experimentally veri- 
fied our theoretical prediction that the diameter of E22 
photoexcited nanotubes initially increases. However, we 
were unable to get a good sample for verifying our re- 
lated prediction that the diameter of E22 photoexcited 
mod 1 nanotubes initially decreases. This will be one of 
the goals of our future studies. 



APPENDIX A: ELECTRONIC STATES 

The tight-binding wavefunctions for the electronic 
states in carbon nanotubes are expanded in terms of the 
symmetry-adapted basis functions 



(Al) 



where r — A, B labels the atoms in the two-atom unit 
cell, Cr{k,^) are expansion coefficients, and |fc,/i, r) are 
symmetry-adapted basis functions. 

The symmetry- adapted basis functions are linear com- 
binations of localized atomic n orbitals 



|/c,^,r) 



(A2) 



where N is the number of two-atom unit cells in the 
system and | J, r) is a localized atomic tt orbital on atom 
r in the two-atom unit cell at Rj. The positions of the 
two-atom unit cells (in unrolled graphene xy coordinates) 
are 



R.i 



where J 



Ji ai + ]2 &2 

-via 



and a2 



/ V3a 
V 2 



(A3) 



(ii,j2) andai = ( - 2 
are the graphene basis vectors. 

The two dimensional wavevector k(/c, ^) appearing 
in the symmetry-adapted basis function expansion in 
Eq. (|A2p is determined by imposing translational and ro- 
tational boundary conditions on the nanotube. Imposing 
the translational boundary condition we have 



k(A:,^) • T = fc T 



(A4) 



where T = ti ai -|-i2 ^2 is the nanotube translational vec- 
tor which is parallel to the tube axis and has the length 
of the translational unit cell. Explicit expressions for ti 
and t2 in terms of the chiral indices n and m can be found 
in Ref. [H and are given by ti = 222±» and ta = 
where du = gcd(2n + m, 2m + n). The rotational bound- 
ary condition is 



k(fc, /i) • Ch = 27r/i 



(A5) 



where Ch = n ai -I- m a2 is the chiral vector in unrolled 
graphene coordinates and = • • • A^i^gx"-'^ angular 
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momentum quantum number that labels the cutting lines 
in the simple zone folding picture,^ 

From Eqs. (jAHIA5|) . we arrive at the symmetry- 
adapted tight-binding wavcfunction 



(A6) 



with the phase factor (/)j(fc, /i) = k(fc, /i) • Rj given by 
7r^((2n + m)ji + (2m -|- n)j2)) 



-t- 



+ nm + 
V3 a k mji — nj2 



y/n? + nm + 



(A7) 



Substituting the symmetry-adapted tight-binding 
wavefunction (|A6[) into the Schrodinger equation, we ob- 
tain, for each value of /x, a 2 x 2 matrix eigenvalue equa- 
tion for the electronic energies Egfj, {k) , and the expansion 
coefficients Cr{s, ^, k), namely 



where s = w, c labels the valence and conduction band 
states. The 2x2 Hamiltonian and overlap matrices are 
given by 

Hry = ^e*'^J'('=^^)(0,r|i/|J',r') (A9) 
J' 



and 



J' 



e"/'j'(fe,f)(o,r|J',r') 



(AlO) 



In the sum over J' ^2)1 in Eqs. (|X9)) and (jAlOp we 

keep only the on-site and third nearest neighbor contri- 
butions since the parameterized matrix elements vanish 
at distances beyond the third nearest neighbor distance 
in graphene. 

Values of (ji, j2) for the first to fourth nearest neigh- 
bors of the A and B atoms in the two-atom unit cell 
J = (0,0) are easy to work out and can be found in 
Table 2 of Ref. M. 



APPENDIX B: PHONON MODES 

Here we derive the block diagonal 6x6 dynamical ma- 
trices and symmetry-adapted atomic displacement vec- 
tors for phonon modes in a carbon nanotube. 

The equilibrium position of the r-th atom in the J- 
th unit cell is R^j = Rj + r^, where ta — 
and Tg — (^^lO). The position of the two-atom unit 
cell in unrolled graphene xy coordinates, Rr-j, is defined 



by Eq. I|A3|) in Appendix [Xj The positions of the near- 
est neighbor atoms are R"j, where a denotes the near- 
est neighbor shell and I labels the atoms within the a- 
th shell. The corresponding instantaneous atomic dis- 
placements are U^j and U"j. In this notation the bond 
stretching potential is given by^i 



E 

rJ.ctl 



pa 
^ s 

2 



(hfj-iV^j^-V^j))' (Bl) 



where is a unit vector parallel to R"j — R^j and F" 
are the bond stretching force constants. 

The nearest neighbor in-plane bond bending potential 
is given by^ 



rJ /(' 



u 



rJ, 



(B2) 



where Tbt is the in-plane bond bending force constant. 
The prime on the sum over nearest neighbors / and V 
indicates that the / = I' term is omitted. 

The out-of-surface bond bending potential is given 



Vrh = 



rJ 



(B3) 



where F^f, is the out-of-plane bond bending force con- 
stant. The unit vector is given by 



QrJ 



E,(R^!,-z (z.R,,j)) 



■R.j))l 



(B4) 



and for large diameter nanotubes is approximately equal 
to the radial outward unit vector at R^j. As pointed 
out in Ref. 'Sl' the form of Qrj in Eq. (|B4[) is required to 
preserve rigid rotational invariance. 

The bond twisting potential is given by^ 



~2~ 



((Ui+U2)-(U3 + U4)) (B5) 



where Ft^, is the bond twisting force constant. In the 
sum represents a bond between atom i and one of 
its nearest neighbors j and f (i.j) is the radially outward 
unit vector at the midpoint of the bond. The atomic 
displacements Ui . . . U4 are the displacements of the four 
atoms attached to the bond with Ui related to U2 
by a C2 rotation about f j) . 

The equations of motion for the carbon atoms are 



M 



rJ 



9t2 



dV 



rJ 



(B6) 



where M is the mass of a carbon atom and V = Vg 
Vbb + Vrb + Vtw is the ionic vibrational potential. 
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Symmetry-adapted atomic displacement vectors, for 
each value of v, can be written as 



i(0j(ij,i^)+cjt) 



(B7) 



where q is the phonon wavevector, v = . . . N-^^^ — 1 
labels the cutting lines, and lo is the phonon frequency. 
The unit vector {q, v) is the phonon mode polarization 
for the r-th atom in the J = (0, 0) two-atom unit cell. 
The phase factor <j)j{q,i') is obtained from Eq. (|A7|) by 
replacing with g, ly. 

The matrix S{9j) is a unitary rotation matrix for ro- 
tations about the nanotube axis (taken to be z), and 
it acts to rotate er{q,v) to the J-th two-atom unit cell. 
The counter-clockwise rotation angle is 9j = ji 6i +j2 O2, 
where 

_ {2n + rnj-K 



and 



(B8) 



(B9) 



+ nm + 



(2m -|- 11)71 
^ -|- nm + m? 



Substituting the atomic displacement vectors (jB7p into 
the equations of motion (|B6|) . we obtain an eigenvalue 
problem for the phonon frequencies and mode polariza- 
tion vectors 

^ Tiryiq, v) er{f3, q, v) = Mu;j^{q) e^(/3, q, v) (BIO) 



where 'Dr,r'{q,v) is a 6 x 6 dynamical matrix for each 
value of V and ujp,j{q) is the phonon dispersion relation 
with (3 — 1 ... 6. Explicit expressions for the symmetry- 
adapted 6x6 dynamical matrix in terms of the atomic 
force constants can be found in Refs. 57 and 58. 

The force constants obtained from the fit shown in 
Fig. [71 are V] = 237.6, V] = 17.47, ^ 2.895, = 
7.166, Tbb = 18.37, E^ = 16.67 and Et„ = 6.609 in units 
of Joules/meter^. 



APPENDIX C: ELECTRON-PHONON 
INTERACTION MATRIX ELEMENTS 



Eq. (jA7p . and er{P, q, are the phonon polarization vec- 
tors in Eq. (|BT01) . In Eq. ([UT]), we choose J" = so that 
S{Oj") — 1 and restrict sums over rJ and r'j' to fourth 
nearest neighbors of Rr".o- 

The deformation potential vectors X{r'J',r",rJ) ap- 
pearing in Eq. (jCip are three-center integrals defined as 



X = J dr ip*,j,{r - Rr'j') Vwc(r - Rr"o) <y5rj(r - R^j) 

(C2) 

where iprji'r — Rrj) is a tt orbital |rJ) localized at R^j 
and Wc(r — Rr"o) is a carbon atom potential centered at 

Rr"0. 

We evaluate the deformation potential vectors A using 
the 2pz atomic wavefunctions and the screened atomic 
potential for carbon in Ref. [l^ obtained from an ab initio 
calculation in graphenci^ The tt orbitals at R^j are 



'^rj(r) = (r • Prj)J2^^ (~^) 



(C3) 



where prj is a unit vector normal to the surface of the 
nanotube at R^ j . Using Eq. (|C3p for the tt orbitals in the 
evaluation of A allows us to take the relative orientation 
of the TT orbitals on different sites into account. Similarly 
the atomic potentials are taken to be 



7^12 (C4) 



Values of vi, ui and r/ are tabulated in Table I of 
Ref. [H. Substituting the expansions (jC3p and (|C4[) into 
Eq. (|C2[) we obtain an expansion for A in terms of three- 
dimensional integrals. In the case of graphene, the three- 
dimensional integrals can be done analyticallyj ^^'^^ For 
the nanotube case we can reduce the three-dimensional 
integrals to one-dimensional integrals that can be done 
numerically. 



The electron-phonon interaction matrix element in 
Eq. ^ is given by 

r" 

C;,{s\p + u,k + q) e-"^3'{k+q.4^+^) 

r'3' 

X C,{s, /i, k) e''^-i(fe.A') A(r'J', r", rJ) (CI) 

rJ 

where Cr{s, /i, k) are the electronic expansion coefficients 
in Eq. (|A8[) . (/)j(fc,/i) is the phase function defined in 
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